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We calculate the quantum radiation power of black holes which are asymptotic to the Einstein- 
de Sitter universe at spatial and null infinities. We consider two limiting mass accretion scenarios, 
no accretion and significant accretion. We find that the radiation power strongly depends on not 
only the asymptotic condition but also the mass accretion scenario. For the no accretion case, 
we consider the Einstein-Straus solution, where a black hole of constant mass resides in the dust 
Friedmann universe. We find negative cosmological correction besides the expected redshift factor. 
This is given in terms of the cubic root of ratio in size of the black hole to the cosmological horizon, so 
that it is currently of order f 0~ 5 (Af/lO 6 M ) 1/3 (t/14 Gyr)~ 1/3 but could have been significant at the 
formation epoch of primordial black holes. Due to the cosmological effects, this black hole has not 
settled down to an equilibrium state. This cosmological correction may be interpreted in an analogy 
with the radiation from a moving mirror in a flat spacetime. For the significant accretion case, we 
consider the Sultana-Dyer solution, where a black hole tends to increase its mass in proportion to 
the cosmological scale factor. In this model, we find that the radiation power is apparently the same 
as the Hawking radiation from the Schwarzschild black hole of which mass is that of the growing 
mass at each moment. Hence, the energy loss rate decreases and tends to vanish as time proceeds. 
Consequently, the energy loss due to evaporation is insignificant compared to huge mass accretion 
onto the black hole. Based on this model, we propose a definition of quasi-equilibrium temperature 
for general conformal stationary black holes. 

PACS numbers: 04.70.Dy, 04.70.Bw 
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I. INTRODUCTION 



Exactly speaking, in our Universe there is no black hole which is asymptotically flat. We call black holes asymptotic 
to an expanding universe cosmological black holes. A particularly important example of such black holes is primordial 
black holes [1] which may have formed in the early universe. They would play a unique role as probes into currently 
unknown physics in various aspects. The key observable phenomenon is the Hawking radiation. See J3] and references 
therein. The effect of cosmological expansion on the dynamics of black holes has been studied 3]. However, the 
cosmological effect on the Hawking radiation has not been seriously investigated yet. It is generally believed that 
the cosmological expansion will not affect the evaporation process if the cosmological horizon is much larger than the 
black hole horizon. Although this is very plausible from a physical point of view, this should be verified from a definite 
argument and it is also important to estimate the possible cosmological corrections. Moreover, this assumption might 
not be the case in some cosmological situations. For example, the black hole horizon may have been of the same 
order as the cosmological horizon immediately after primordial black holes formed and/or if matters around black 
hole continued accreting so rapidly that self-similar growth of black hole horizon might be possible [H . 

The evaporation of cosmological black holes is also important from a point of view of thermodynamics. The 
black hole thermodynamics [4| offers a unique possibility of understanding the theory of gravity through the laws of 
thermodynamics. It is argued that a black hole is described as an object in thermal equilibrium (black body) with 
temperature Th = k/2tt [a, @, where K is the surface gravity of the event horizon, and evaporates by radiating its 
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mass energy according to the Stefan-Boltzmann law [7| ■ This argument has been so far established only for black holes 
which are asymptotic to flat, dS and AdS spacetimes at spatial and null infinities. See [4( for asymptotically flat case, 
and [|[ for asymptotically dS and AdS cases. It might be reasonable that black hole thermodynamics requires the 
asymptotically static nature of spacetimes, because an equilibrium state corresponds to a static condition. However, 
if black holes can be regarded as thermodynamic systems even for nonstationary cases, it is natural to ask what kind 
of nonequilibrium states correspond to dynamical black holes. An interesting example is cosmological black holes. 
The spectrum of its Hawking radiation has already been studied for the case of no mass accretion [9|. However, the 
radiation power or luminosity has not been explicitly calculated yet. 

In the present article, we calculate the power of the Hawking radiation through the quantum expectation value of 
stress-energy tensor in an expanding universe based on the two limiting scenarios about mass accretion. To do this, 
we use two interesting models for cosmological black holes. One is the Einstein-Straus black hole and the other is the 
Sultana-Dyer black hole. The former has no mass accretion, while the latter has significant mass accretion. Our first 
model was raised by Einstein and Straus [10j. This is the exact solution of the Einstein equation with timelike dust, 
which is obtained by pasting a Schwarzschild spacetime with a dust Friedmann universe on a timelike hypersurface. 
This kind of models are often termed as of "Swiss-cheese" type. The existence of a black hole event horizon in the 
Einstein-Straus solution is guaranteed by construction and the radius of the event horizon is constant. All energy 
conditions are of course satisfied in this spacetime. Our second model was given by Sultana and Dyer [l2|, This is 
obtained by a conformal transformation operated on the Schwarzschild spacetime, whose conformal factor is carefully 
chosen so that the metric is the exact solution with the combination of timelike and null dusts and the spacetime 
is asymptotic to the Einstein-de Sitter (or flat dust Friedmann) universe at spatial and null infinities. Since this 
transformation does not affect the causal structure, the existence of an event horizon is guaranteed, although this 
spacetime has some trouble with energy conditions. This spacetime has a feature that the physical radius of the event 
horizon increases due to accretion, and it approaches infinity as time proceeds although its growth rate tends to be 
much slower than the growth rate of the cosmological horizon. 

This paper is organised as follows. Sections |TT] and IIIII are devoted for the calculation of Hawking radiation 
from Einstein-Straus and Sultana-Dyer black holes, respectively. Summary and discussions are given in section 
ITVl Throughout this paper, we use the Planck units, c = ^ = G = fcs = l. 



II. HAWKING RADIATION FROM THE EINSTEIN-STRAUS BLACK HOLE 



A. The Einstein-Straus black hole and its formation 



where 



The Einstein-Straus black hole is constructed by pasting the Schwarzschild and the Friedmann solutions at a 
spherically symmetric timelike hypersurface £ (see [1G | or Appendix A in Q). 
The Schwarzschild metric is given by 

ds 2 BB = -C(R)dT 2 + -^ + R 2 dn 2 , (2.1) 
2M 

C(R):=1- — , (2.2) 

M is the mass of this black hole, dQ is the line element of a unit two dimensional sphere, and T and R are the time 
coordinate and areal radius in the Schwarzschild coordinates, respectively. We can get the double null form of the 
metric as 

rfs BH = -C{R)dUdV + R 2 dn 2 , (2.3) 

where 

«■ : = mr (24a> 

U:=T-R*, V:=T + R*. (2.4b) 



The Friedmann metric is given by 



H = ~<H- "f'r ( j^+r 2 d^ j . (2.:,) 
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where a is the scale factor, k is the spatial curvature, r is the comoving radius and t is the proper time of the comoving 
observer or the cosmological time. Hereafter we assume the open or flat Friedmann metric (k = or —1) in order to 
guarantee the existence of future null infinity. The double null form of the metric is given by 

ds% = -a 2 du dv + R 2 dn 2 , (2.6) 

where 

dt dr . 

dv ■= — , d X := , (2.7a) 

a VI- kr 2 

u := r/—x, v:—rj + x- (2.7b) 

The junction surface E is given by a constant comoving radius r = and located at 

R = a(t)rz- (2.8) 

The Israel junction condition [Tl| with no singular hypersurface reduces the continuity of the first and second funda- 
mental forms on the junction surface E. Then, we obtain two independent equations, 



where a := da(t)/dt and 



dT(t) Jl - kr 2 
a J a 1 r£ a A 



C| s = l- — . (2.10) 



Equation (|2.9a|) gives the relation between the time coordinates T and t. Equation (|2.9b[) determines the time evolution 
of the scale factor a(t) and identical with a dust Friedmann equation. Thus, it is required that the Friedmann metric 
is that of the dust Friedmann universe. By comparing eq. (|2.9bl) with the dust Friedmann equation, 



\ 2 

a \ k 8ir p„ 



a J a 2 3 a 3 
where p* is a constant, we find a relation between p» and M , 



(2.11) 



M = -7rr!/9*. (2.12) 

The Einstein-Straus black hole describes a cosmological black hole spacetime with no mass accretion. 

To calculate the power of Hawking radiation from the Einstein-Straus black hole, we need to specify how it has 
formed in the expanding universe. Here we assume that an overdense region of which comoving radius is rs and 
mass M begins to contract at the moment t = is and collapse to form a black hole. The Penrose diagram of this 
gravitational collapse is shown in fig. [1] In this figure, H, I, B, S denote the event horizon, the null infinity, the 
spacelike hypersurface t = *b for r < and the surface of the collapsing dust ball, respectively. We can assume that 
the collapsing region surrounded by B, S and the regular centre is described by some regular dynamical metric, 

ds 2 col = A(t, A) (-dr 2 + d\ 2 ) + R(t, X) 2 dil 2 , (2.13) 

where r and A are respectively appropriate temporal and radial coordinates. We can set that A = corresponds to 
the regular centre. The double null form of the metric is given by 

ds 2 col = -Adad(3 + R 2 dQ 2 , (2.14) 

in the collapsing region, where 

a:=T-X, f3:=r + X. (2.15) 
As assumed above, the starting surface B of collapse is given by a spacelike hypersurface in the expanding region, 

t = t B {r) and < r < r s . (2.16) 
This surface is also described in the collapsing region by 

T = r B (r) and A = A B (r) , (2.17) 
using the coordinates in the collapsing region. 
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FIG. 1: The Einstein-Straus black hole formed in an expanding universe. The zig-zag lines denote spacetime singularities. H, 
I, B, S and E denote the black hole event horizon, the null infinity, the starting spacelike surface of gravitational collapse, 
the surface of collapsing dust ball and the junction surface between the Schwarzschild and Friedmann spacetimes, respectively. 
The empty region surrounded by 5", E and the black hole singularity is described by the Schwarzschild solution. The collapsing 
region surrounded by B, S and the regular centre is described by some regular dynamical metric. The expanding region 
surrounded by B, E, I, the regular centre and the big-bang singularity is described by the Friedmann solution with dust. 
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FIG. 2: The Penrose diagrams of the cosmological black hole with no matter accretion. Upper panel: The initial surface and 
the comoving observer at r b s > rs are shown. The null geodesic 7 passes near the event horizon. A late-time observer detects 
a scalar wave which has propagated along 7. Lower panel: The Penrose diagram of gravitational collapse extended to the 
negative radius region. This extension makes it easier to calculate (T^). 



B. Redshift and Hawking radiation 



We introduce a matter field (f> which describes quantum radiation from a black hole. For simplicity, let be a 
massless scalar field with minimal coupling, which satisfies the Klein-Gordon equation, □</> = 0. In manipulating 
quantum field theory in curved spacetimes, especially on black hole spacetimes, we need to estimate the redshift. 

The wave mode of 4> propagates along a null geodesic 7 passing near the event horizon. As shown in the upper 
panel in fig. [21 this mode is ingoing at the initial surface and becomes outgoing after passing through the centre. The 
function v = G(u) relates the ingoing null coordinate v of 7 at the initial surface and the outgoing null coordinate u 
of 7 at late times. Here note that, although both the initial surface and the spacetime region where a distant observer 
is are given by the Friedmann solution, we need to distinguish the null coordinates by (u, v) at the initial surface and 
by (u, v) for the observer because their values are different from each other in this construction. The function G(u) 
is obtained by the junction of null coordinates at the intersections of 7 with the surfaces B, S and E. 
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The junction at B gives a relation between v and (3 as v = v((3). Since we can assume that the gravitational collapse 
begins with sufficient smoothness, v = v{(3) is C 1 for the relevant ingoing null rays, i.e., 



v — vi (3 + vq , 



(2.18) 



where vq and vi(> 0) are constants [22|]. The reflection of 7 at the regular centre is given by a simple replacement of 
(3 by a as v = v(a). This expresses the redshift along 7 from the initial surface to the collapsing region. 

The junction at S gives a relation between a and U as a = a(£/). This junction is completely the same as one gets 
for an asymptotically flat black hole (see Q and §8.1 in [(|). Hence, we do not show the details of the calculation but 
only quote the result: 



ct\ exp (—k U) + «o • 



where ap and «i(> 0) are constants and 



ldC 
2dR 



R=2M 



1 

AM' 



(2.19) 



(2.20) 



This expresses the redshift from the collapsing region to the empty region. 

The junction on E gives a relation between U and u as U — U(u). To obtain this, we need to consider the junction 
of metrics, ds| H | s = dsp| s . Since X is the timelike hypersurface given by r = in the comoving coordinates, the 
null coordinates on S are regarded as functions of t. Furthermore, from the relation (|2. 7a|l . these null coordinates are 
instead regarded as functions of rj. Then the junction of the null coordinates on £ can be discussed using the partial 
derivative with respect to rj, i.e., the junction of metrics gives the following relation on S, 



(2.21) 



On the other hand, from eqs. (|2.4[) and (|2.8p . we get on T, 



V ,7) U ,7} — Z 



(2.22) 



where the prime ' denotes the argument differential, a' := da(rj)/dr]. Therefore, eq. (|2.21[) becomes a quadratic 
equation for U :V on S. We get the positive root and obtain dll/du as 



dU 
da 



U, 



i 2 C 



(2.23) 



This expresses the redshift along 7 from the empty region to a comoving observer in the expanding region. 
Thus, we obtain the function G(u) by combining three functions and the reflection at the centre. 

v = G(u) — v[ exp [—k U (u)] + Vq , 



(2.24) 



where v' and w^(> 0) are constants. Furthermore, for later convenience, we will extend the background to include 
the negative radial coordinate region. The extended Penrose diagram is shown in the lower panel in fig. [5] Then we 
do not need to consider the reflection of 7 at the centre. By this virtual extension, eq. (|2.18[) becomes u — u\a + u' Q , 
where u' and ui(> 0) are constants. Hence, the redshift in the extended background spacetime is given by the same 
function G(u) obtained above with replacing v by u in the left-hand side, 



= G(u) = u\ exp [— k U (u)] + uq , 



(2.25) 



where u is a constant. Note that this function (|2.25p also applies to the dimensionally reduced spacetime introduced 
below. 



C. Hawking radiation from cosmological black hole with no accretion 



The quantum expectation value must be renormalised. The regularisation technique in two dimensions has been 
well established. Appendix summarises the calculation of the vacuum expectation value of stress-energy tensor 
(vac I I vac) and its renormalised value {T^ v ), where |vac) is an appropriate initial vacuum state. For simplicity, 
we here reduce the gravitational collapse spacetime described by fig. |5]to a two dimensional one by cutting out 
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the two dimensional angular part from metrics (|2.ip . (|2.5[) and (|2.13[) . It has been known that two dimensional 
Schwarzschild black hole gives a qualitatively correct power of Hawking radiation in four dimensions, since the so- 
called grey body factor in the Hawking radiation disappears due to the absence of curvature scattering of matter fields 
in two dimensions. The curvature scattering of matter fields also does not occur in two dimensional Einstein-Straus 
spacetime. We can expect that two dimensional Einstein-Straus black hole gives a qualitatively correct radiation 
power. 

The thermal radiation in asymptotically flat black hole spacetimes has been obtained under the following three 
procedures [H, [|| ; to neglect the curvature scattering, to define an initial vacuum state on a spacelike hypersurface 
before the black hole formation, and to observe particles at sufficiently late times. The first is automatically done if 
we work in two dimensions. The second implies that, in the Hciscnbcrg picture, the quantum state which has been 
vacuum initially is no longer vacuum after the gravitational collapse. The third implies that a wave mode detected by 
a distant observer should pass the neighborhood of the event horizon and hence it has been strongly redshifted before 
it is observed. This means that it has been of very high frequency near the event horizon, and the geometrical optics 
approximation is valid, being consistent with neglecting the curvature scattering even for four dimensional case. 

To calculate the Hawking radiation from the Einstein-Straus black hole, we adopt the vacuum state |vac) associated 
with a comoving observer at the initial surface as a physical initial vacuum state and calculate {T^ v ) for a distant 
comoving observer at sufficiently late times. To be precise, |vac) is defined by the quantisation of (f> using the normal 
modes obtained in the coordinate system of eq. (|2.7b[) on the initial surface, and the components of (T^ v ) is calculated 
in the same coordinates. Moreover, in two dimensions there is no genuine cosmological particle creation for a massless 
scalar field (see §3.4 in [6j for example), and hence (T^) expresses purely the Hawking radiation from cosmological 
black holes. 

First, we calculate (T^ p )( in ' at the initial surface using eqs. (|A.4[1 and (|A.5[) given in Appendix. The metric suitable 
for this purpose is given by ds 2 ? = —a 2 dudv, for which 

2[(a') 2 - aa"] 

U = — — — - (2.26) 



in eq. (|A.4[) and D — a 2 in eq. (|A.5[) . Therefore we obtain 



/ T __v(in) _ , T v(in) _ 2 ( a ') -a a" 

48n or 



(2.27a) 



(T _ \(in) {a') 2 -a a" . 

The observed power for a comoving observer on the initial surface is proportional to fj-x component of (T^)^. By 
coordinate transformation from (u, v) to (fj, x) > we nn d 

(T m )^ =0. (2.28) 

This indicates that no energy flux is observed at the initial surface. 

Next we calculate (7) ll /) ( - obs ' > measured by a comoving observer at late times. The metric suitable for this purpose is 
ds 2 ? = —a 2 dudv. This gives the same 1Z as for (Tp,p)( ia \ In order to find the function D in eq. (|A~5l) . null coordinates 
(it, v) of a late-time observer should be expressed in terms of null coordinates (u, v) with which the initial vacuum 
state | vac) is defined. The relation between (u, v) and (u, v) reflects the time evolution of the background spacetime 
in between and is given by the rcdshift along outgoing and ingoing null geodesies which connect the late-time observer 
and the initial surface. For the extended background spacetime (lower panel in fig. [5]), the relation between u and 
u is given by the redshift (|2.25[) of outgoing null geodesic 7. The relation between v and v is a simple one v = v, 
because the relevant ingoing null geodesic lies in the expanding region without passing any surfaces B, S and X during 
propagating from the initial surface to the late-time observer. Consequently, D which gives the metric at late times 
as dsp — —Ddudv is obtained as 

„ o du a 2 , 
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where the coordinate transformation (|2.7b[) is used and G' :— dG(u)/du. This gives 

a 2 G" 



a a 

qJ2 - q,3 

(a a')' 

2G' 3 
a a' 
"G 7 " ' 
(aa'Y 

2G' 



G' 4 



i 2 G"' 
G' 4 



2r<ii2 



a 2 G 



G' 5 



where a jfi = (du / 'du)(drj / 'du)a' = a'/2G' and a = (dv / dv)(dr] / dv)a! — a'/2 are used. Hence eq. (|A.4[) gives 

2 



/rrn \(obs) _ ' 



(TV 



i (obs 



\ (obs) 



1 


3 / G" N 


24^ 




2(a 






487T0 2 


(a') 2 


-a a" 


487ra 2 ' 



G'" 
"G 7 



(2.30a) 
(2.30b) 
(2.30c) 
(2.30d) 

(2.31a) 

(2.31b) 
(2.31c) 



The observed power P b s is given by the tetrad component (T^)( obs \ We obtain P b s from the above calculations, 



/rp(x)\(obs) IT \(°bs) 



24tt og 



3 (&[_ 
G' 



G'" 



(2.32) 



where ciq := a(?7o) an< i Vo is the conformal time at the moment of observation. By substituting the expression (|2.25|) 
for G(u) into P bs, 



Pobs — 



EL 

48-7T V a 



24tt al 



El 
W 



U'" 

IF 



(2.33) 



where U' := dll/du is given by eq. (|2.23[) . 

By comparing eq. <|2.33|) with eq. (|2.28|) . it is obvious that the quantum creation of energy flow occurs due to 
the forming black hole. Here, recall that the power of the Hawking radiation Ph(2D) m an asymptotically flat two 
dimensional black hole is 



Pi 



H(2D) 



1 

2^ 



<ko- 



exp(27rw / k) — 1 

\2 • 



K 

48^ 



(2.34) 



Comparing P h s with Ph(2D)? we find that the factor (U'/ao) in the first term and the whole of the second term in 
eq. (|2.33p are the effects of cosmological expansion. 

Here we should recall that our calculation is performed on a two dimensional background spacetime. That is, in 
calculating Ph(2D) m ec l- (|2.34[) , the state density N/2tt at energy level u> is appropriate to two dimensional case (one 
spatial dimension), where N is the effective degrees of freedom and N = 1 for a scalar field. Therefore the numerical 
factor in eq. (|2.33p will be valid only for two dimensional case. However, we expect that eq. (j2.33[) qualitatively correct 
even for four dimensional case if we neglect the curvature scattering and the cosmological particle creation. 



D. Application to two dimensional Einstein-Straus black hole 

Here we apply eq. (|2.33[1 to our collapse model shown in fig. [2] In the following, we assume k — for simplicity. 
Then, from eq. (|2.9b[) and the relation adrj = dt, the scale factor becomes 




(2.35) 
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where r)i n = 3t; n , and i; n and r)i n are respectively the cosmological and conformal times at the initial surface. We 
normalise the scale factor at the initial surface. Furthermore, the Friedmann equation (|2.11[) relates t- m and rj ln with 

P*, 



P* 



2tt< 



Then eq. (j2~T2l) gives 

M = — . 

Substituting eq. (|2.35[) into the right-hand side of eq. I|2.23jl . we obtain 

a 



U' 



F 



where 



F := 1 



2r s 
»7 



where eqs. f|2.37[) and r?i n = 3ii„ are used. Hence, substituting eq. ()2.38[) into eq. (|2.33[> . we obtain 



obs 



a rct 

a 



2 2 



1 



48^ F(r) Te 



Qrct 

a 



247ra r 2 ct 



3 (a! 



a 
a 



~^~F 



El - 1 (EL 
~f " 2 If" 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



(2.40) 



where Q le t denotes the evaluation of Q at r\ = r] le t :— rjo — (r b s — ru) when the ray 7 intersects E. This P bs is 
regarded as a function of the cosmological time to of the observer by using eq. (|2.35|) , which gives 



, ^in 3 ^in 

tret = T7T ^rct = ~T 



V 



'/in ,1/3 / x 
"1/3*0 -(^obs-^s) 



f 



,1/3 



G 



M\ 1/3 /rob £ _ 1 



(2.41) 



where eq. (12.371) is used in the last equality. 

Furthermore the observed power P ^ s can be expressed in a more convenient form. Using eq. (|2.37[) and the Hubble 
parameter 



H 



we can get 



2rs 

Vret 



(2MH ret ) 



2 

r/a ' 

1/3 _ ,1/3 



(2.42) 



(2.43) 



where e := 2MH rct is the ratio of the black hole horizon radius to the Hubble horizon radius when 7 intersects E. 
Furthermore, we observe the cosmological redshift z of the photon emitted from E, 



1 + z 



&rct 



(2.44) 



This z will be regarded as the redshift of the host galaxy of the black hole, and the ratio e can be expressed as 
e = 2MHq{\ + z) 3 / 2 , where Hq is the present Hubble parameter. It is very natural that the cosmological correction 
is given in terms of the ratio e. Using this ratio, we can express P D b s in eq. (|2.40[) simply as 



-fobs — 



48tt (1 + zf 



(1 



J./3N-2 



A/3 



1 + e 1 / 3 8 Kl + e 1 / 



cl/3 



(2.45) 



where k = 1/(4M). 
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Note that the observed power P b s is not intrinsic but cosmologically redshifted. The intrinsic power Pes(2D) is 
then given by 



ES(2D) 



(1 



(1 



,1/3 



r 2 + < 



,1/3 



,1/3 



,1/3 



,1/3 



(2.46) 



The evaporation should be described by this intrinsic power. Up to 0(e 1 ^ 3 ), we get -Pes(2D) — ^h(2D) (1 — 26 1 / 3 ). 
This implies that the intrinsic power is suppressed by the cosmological expansion. The physical interpretation of this 
effect is proposed in section IIVI 

We are interested in two distinct limits from a physical point of view. In the first, the event horizon is much smaller 
than the Hubble horizon at present. This corresponds to the limit e — > with keeping z constant, and we obtain 

p r, 17-> r, //I . V> 1 .1 - -1 

*ES(2D) _ 

evolution. 



* ^H(2D) 

i.e., r] 



and P bs — > -Ph(2D)/(1 + z ) 2 ■ In the second, we consider a very late phase of the cosmological 
-> oo. This corresponds to the limit e — > and z — > simultaneously as seen from eqs. (|2.42|l and 



(|2.44|) . Then we obtain P ES (2D) -> -Ph(2D) and P Q b s -> Ph(2D)- 

In black hole thermodynamics the Schwarzschild black hole is regarded as in thermal equilibrium, and the 
temperature Th is assigned to the black hole (zeroth law). This temperature is given by Th = k/2tt — 1/(8ttM) which 
satisfies the Stefan-Boltzmann law in two dimensions Pr(2D) — (7r/12) T§ as seen from eq. (|2.34l) . Then one might 
also want to assign the temperature Th = l/(8irM) also to the Einstein-Straus black hole. However, the radiation 
power Pes(2D) deviates from the Stefan-Boltzmann law due to the correction term of 0(e 1//3 ). This suggests that the 
Einstein-Straus black hole deviates from thermal equilibrium in a finite cosmological time. Only in the limit 770 — > 00, 
this black hole settles down to thermal equilibrium. 



E. Evaporation of the Einstein-Straus black hole in four dimensions 

The power Ph(2D) = n 2 /A8tt is obtained for asymptotically flat two dimensional black holes. Therefore, we simply 
replace the factor k 2 /487t in eq. (|2.46[) by the four dimensional counterpart Ph(4D)- This Ph(4D) is given by the 
Stefan-Boltzmann law in four dimensions, 



P 



H(4D) 



aT^An = 



N 



30720ttAP 



(2.47) 



where An = 4ir(2M) 2 is the spatial area of the event horizon and a — Ntt 2 /120 is the Stefan-Boltzmann constant for 
the massless matter field with the effective degrees of freedom N. Here TV is given by 
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N := rib H — n i ■ 



(2.48) 



where rib and rif are the numbers of helicities of massless bosonic and fermionic fields, respectively, and the factor 
7/8 comes from the difference of statistics of fermions from bosons (see for example [7J for derivation). Then it is 
appropriate to estimate the order of iV by the standard particles (inner states of quarks, leptons and gauge particles 
of four fundamental interactions), N ~ 100 if the black hole temperature is lower than ~ 1 TeV. Next we consider 
about the correction terms in square brackets in eq. (|2.46|) . These terms come from the factors U', U" and U'" in 
eq. (I2.33[) . Here recall that the function G(u) in eq. (|2.25[) is valid even for four dimensions. Hence we can expect 
that the same correction terms appear as well for four dimensional case. From the above consideration, we expect 
that the four dimensional intrinsic power Pes(4D) is given by 



P 



N 



es(4D) 3o 7 207rM 2 



(1 



,1/3^-2 



,1/3 



)- 2 + 8eUl + 



1 



1 + e 1 / 3 8V1 + C 1 / 



,1/3 



(2.49) 



Finally we estimate the evaporation time of the Einstein-Straus black hole £es- Equations (|2 .35[) and (|2.43[) give 
e = (4M)/(3i re t)- Hence, equating Pes(4D) to —dM/dt re t in the left-hand side of eq. (|2.49|) . we can regard it as the 
evolution equation of mass M as a function of the cosmological time t. Up to the first correction term of order (^(e 1 / 3 ), 
eq. (|2.49[) gives the semi-classical evolution equation of M(t Ie t) as 



dM 
~dT 



N 



307207rAf 2 



2 (") 



1/3- 



(2.50) 
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where we denote t rc t as t, representing the cosmological time of the evaporating black hole. Since the correction is 
negative, the emission is suppressed and the life time is prolonged. Assuming that the correction term is small, we 
get the order estimate for the deviation of £es from the evaporation time of the Schwarzschild black hole in as 

— - 1 = 0(A/- 2 / 3 ) , (2.51) 

or 

tES -tn = 0(M J / 3 ) , (2.52) 
where % is given by neglecting the correction term in eq. (|2.50p , 

30720tt , r , 

tn =s M ■ (2.53) 

We can see from eq. (|2 . 5 1 1) that as the initial mass is larger, the evaporation time is better estimated by £h- This 
is reasonable since the Hubble parameter of the Einstein-de Sitter universe becomes small as time proceeds and the 
cosmological effect on the evaporation becomes negligible. It should be noted that the cosmological correction on the 
evaporation time is relatively small even for a primordial black hole even if it was as large as the particle horizon 
unless its mass is of order the Planck mass. On the other hand, the deviation (tgs — *h) itself can be very large if the 
black hole is very massive. 

III. HAWKING RADIATION FROM THE SULTANA-DYER BLACK HOLE 

A. The Sultana-Dyer black hole 

The Sultana-Dyer black hole is obtained by the conformal transformation of the Schwarzschild black hole [T3|. Its 
metric is given by 



-drj 2 + dr 2 + r 2 dfl 2 + —(drj + dr) 2 
r 



(3.1) 



where M is a positive constant, a(jf) = {rj/rj^,) 2 and 77* is a constant. This spacetime is asymptotic to the Einstein- 
de Sitter universe as r — > 00. Here we consider the following coordinate transformation, 

r ? :=< + 2Mln(^--l) . (3.2) 



This transforms the metric (|3.1[) to the conformal Schwarzschild one, 



dsl D = a(t,r) 2 



r I \ r 



(3.3) 



r = 2M remains an event horizon because the conformal transformation preserves the causal structure. The Penrose 
diagram of this spacetime is shown in fig. [3] There are curvature singularities at r\ = and r = 0. The singularity at 
77 = is spacelike for r > 2M, timelike for r < 2M and null for r = 2M. The central singularity at r — is spacelike 
and surrounded by the event horizon. Hereafter we consider the spacetime given by regions I and II shown separately 
in the right panel of fig. [3] 

This spacetime is conformally static, since there exists a conformal Killing vector £ = dt which is the Killing vector 
on the Schwarzschild spacetime and satisfies the following relation 

£f 9nu = (£j In a 2 ) g^y , (3.4) 

where dt = d n due to the coordinate transformation (|3 - 2[) and £^lna 2 = 4/77. The hypersurface at r = 2M is the 
conformal Killing horizon which is a null hypersurface where £ becomes null. This coincides with the event horizon 
of the Sultana-Dyer black hole. The Misner-Sharp mass m at an arbitrary spacetime point is given by 

m(7y, r) = Ma - 2Mra' + L^LL h + (3.5) 

2a \ r J 
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FIG. 3: The global structure of the Sultana-Dyer solution. The zig-zag lines denote spacetime singularities, while the solid 
lines denote null infinity. So, the three pairs of regions, [I and II], [III, IV, V and VI] and [VII and VIII], respectively describe 
universes disconnected from each other. Geodesies can reach the line of r\ — -co and r = 2M within finite affine length. The 
spacetime given by Regions I and II is asymptotic to an expanding Friedmann region and represents a cosmological black hole. 



where a' := daj dr\ — 2r\jrfc. Then the Misner-Sharp mass at the event horizon is 

„, / 8M 32M 2 \ , o „, 

m EH := Ma 1 + — — . (3.6) 

V V T J 

This means that the mass of the event horizon tends to increase in proportion to the scale factor a cx rf as r/ — > oo. 
Substituting the metric (|3.ip into the Einstein equation, we get 

T ^l D) = Pm u v + p r fc M k v , (3.7) 

where is a unit timelike vector, k^ is a null vector normalised by k^u^ — — 1, p m is the density of timclikc dust 
and p-c is the density of null dust. The vectors and k^ are given in (77, r, 0, ip) coordinates as 



' r rj 2 y/r 2 + 2M(r - rj) r r] 2 ^r 2 + 2M(r - 77) 



The densities of dusts are given by 



4 12[r 2 + 2M(r-r 1 )] 

p™ = ^72^6 ' ( 3 - 9a ) 

4 AM[4r 2 + 3M(2r- V )] 
Pl V * 8tt r 2 r, 5 [r 2 + 2M(r - rj)] ' [ ' 

These densities imply that energy conditions are satisfied only when r\ < r(r + 2M)/2M where p m > and p T > 0. 
Furthermore the velocities of dusts (|3 . 8[) denote that the null dusts in the region i] < r(r + 2M)/2M fall towards the 
black hole, and also that the timelike dusts in the region 77 < r(r + 2M)/2M and 77 < 2r do so. We can describe a 
physical picture that the accretion of timelike and null dusts increases the mass of black hole as shown in eq. Q3.6P . 
For 77 > r(r + 2M)/2M, the source matter fields of the Einstein equation get unphysical. However, the Sultana-Dyer 
metric is featured with the global structure of a cosmological black hole as seen in Fig.[3]and also with the conformally 
static nature which makes the physical interpretation and calculation of quantum stress-energy tensor most tractable. 
Hence, we adopt it as a not only workable but also physically interesting model for a cosmological black hole with 
significant mass accretion. 
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One may consider the future outer trapping horizon as a local definition of a black hole horizon [14| . The trapping 
horizon is given by 2m(j],r) — R(ri), where R — ra is an areal radius. Thus the trapping horizon in this spacetime is 
obtained by the following algebraic equation: 

2M 8M 4r 2 / 2M 



r rj r\ 2 \ r 



1 + . (3.10) 



This has two roots, r — r\ and r = r 2 (r\ < r^), where 



w , -V+ V V + 24Afr7+16A7 2 
n := —M H , (3.11a) 

r 2 ■= \ (3.11b) 

Regions < r < n, T\ < r < r% and r 2 < r are future trapped, untrapped and past trapped, respectively. 



B. Power of Hawking radiation from the Sultana-Dyer black hole 

We introduce a matter field which describes quantum radiation from the black hole. Let </> be a massless scalar 
field with conformal coupling, satisfying (□ — lZ/6)4> = 0, where 1Z is the Ricci scalar, for which we will see that 
the renormalised stress-energy tensor (T M „) can be expressed in terms of (T MJ/ )s c h for the Schwarzschild spacetime. 
Appendix summarises how to derive (Ijj„). We use eq. (|A.15[) given in Appendix. For the Sultana-Dyer spacetime, the 
metric and the tensors X^ and Y" M „ in eq. (|A.15[) are all for the Schwarzschild spacetime, and (T^ v ) is (T , A , y )sch- 
Obviously, A M „ = = since R^ v = 0. Then eq. (|A.15j) becomes 



where A M „ and Y^ v are obtained by substituting the metric (|3.ip into eqs. (|A. 16|) . The first term in (T^) expresses 
purely the Hawking radiation from black hole, and the second term includes cosmological particle creation. Hereafter 
we consider an observer at r = const, i.e., we calculate (T^) in (77, r, 8, tp) coordinates of eq. (|3.1[) . which implies that 
the vacuum state |vac) in is defined with respect to the mode function of cj> in the same coordinates. 

For the observer distant from the black hole, the observed energy flux F b s is given by 

F ohs := <T«> = -±(T vr ) Sch + 1 (l Xvr _ , (3 . 13 ) 

where (T^) is a tetrad component evaluated in the Sultana-Dyer spacetime, and eq. (|3 . 1 2[) and T^) = —a~ 2 T nr at 
the distant region are used in the second equality. The distant observer is comoving with the timelike dust as seen in 
eq. (|3.8p . From eqs. (|A.16|) . the second term in the right-hand side of eq. (|3.13[) becomes 

2880^2 (e x,,r " Yw 

M 3 rj* M 2 T]\ [-9M rf + (21M + r,)^ + \2Mr)l] 



30tt 2 rj 9 r 6 90?r 2 77 14 r 5 

M 2 T)\ [~%%Mrf + (28M - r/)^ 4 + 160 Mr)*] 
120tt 2 r/ 15 r 4 

(3 14) 

Mi/J [2M(65M - 37?7)77 8 + (-240A/ 2 - 64Afry + Ht; 2 )?? 4 ?? 4 - 120M(5M - 77)77®] v ' ' 

180tt 2 r? 16 r 3 

M rj 4 [-4M(65M - 27rj)rf + 5(96M + llr?)*,^ 4 + I2OOM778] 



360?r 2 77 16 r 2 
Afr; 4 [l3r; 8 - 2Ar) 4 r) 4 - 6O77S] 
727r 2 rj 16 r 

This falls off very rapidly for a distant observer, and hence we get -Fobs = — a 4 (T vr )sch for the distant observer. Here 
the coordinate transformation (|3.2p gives (T^ r )s c h = (Ttr)sch + (1 — r /2M) _1 (T t t)sch- Hence we obtain 



obs 



1 , 1 ,_,(*) 



,4 



(T tr ) S ch = -4 (T{;>) Sch , (3.15) 
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where (T^j)sch is a tetrad component evaluated in the Schwarzschild spacetime. This F b s is the flux (energy flow 
per unit time and unit area) detected by the distant comoving observer. 

It should be pointed out that the flux P bs contains only the Hawking radiation from black hole but no cosmological 
particle creation, since F bs is proportional only to (T tr )sch as shown in eq. (|3 . 1 5|) . This does not imply the absence 
of cosmological particle creation at our distant observer. The energy density a of the quantum field <fi in the distant 
region indicates the cosmological particle creation at the distant observer. We get from eqs. (|3.12[) and (|A.16|) . 



: - ( T (v)) - a 2 ( T vv) 



--At, 



7777/Sch 



100?7 ; 
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(3.16) 



240ttV 



for the distant observer. The second term does not include M and expresses purely the cosmological particle creation 
in the distant region. However this raises no energy flow as shown in eq. (|3.15[) . 

Hence F b s in eq. (|3 . 15|) is the flux of the Hawking radiation from Sultana-Dyer black hole. Then the intrinsic 
power of the Hawking radiation Pgo should be given by 



P SD := 4-7r(ra) 2 F obs = -~ 



47rr 2 (T ( ( ) ! ) ) ) Sch 



(3.17) 



where the right-hand side should be evaluated for the distant observer. Here note that the factor 47rr 2 (T^j)s c h is 
the observed power of the Hawking radiation in the Schwarzschild spacetime, and it should equal Pr-(4D) given in 
eq. (|2.47[) if the geometrical optics approximation is valid. Therefore, under this approximation, we get 



1 Nk 2 



a 



2 1920tt ' 



(3.18) 



where n = I /{AM) and N is given by eq. (|2.48|) . 

The geometrical optics approximation gets very good for late times, and the mass toeh(??) in eq. (|3.6|) becomes 
wieh — * Ma as rj — > oo. Hence, by comparing our result Psd with the Schwarzschild one Ph(4D)> h is suggested that 
the effective temperature T c ff of the Sultana-Dyer black hole at late times is given by 



T. 



1 



off 



in Ma 2tt a 



(3.19) 



So, both the intrinsic power and temperature of the radiation from the Sultana-Dyer black hole arc the same as those 
of the Hawking radiation from the Schwarzschild black hole of which mass is the momentary mass of the growing event 
horizon. This temperature T c g and the intrinsic power Psd decrease as time proceeds. This result is reasonable since 
the Sultana-Dyer black hole describes significant mass accretion. This black hole does not lose but gain mass due to 
the accretion of timelike and null dusts, and can be regarded as an object in quasi-equilibrium with temperature T e ff . 



C. Conformal dynamics at infinity and temperature of Sultana-Dyer black hole 

A stationary spacetime is defined by a timelikc Killing vector £ satisfying Cq g^ v = 0. The Killing horizon is a null 
hypersurface where £ becomes null, and the surface gravity n is defined by ('"VaC = at the Killing horizon. 
The value of k changes according to the normalisation of C, by definition. For asymptotically flat stationary black hole 
spacetimes, the Killing horizon coincides with the event horizon. The surface gravity k of the stationary black hole 
is constant everywhere on the event horizon. This is the zeroth law of black hole thermodynamics [J, Q which states 
that a unique temperature can be assigned to the stationary black hole. Then the thermal spectrum of Hawking 
radiation from the stationary black hole, which is a quantum phenomenon, determines the value of the temperature 
to be k/2it under the normalisation of Killing vector as C m Cm ~* — 1 a ^ nuu an d spatial infinities [5|. 

Several generalisations of the zeroth law have already been discussed for general conformal stationary black hole 
spacetimes whose metric g^ v is given by g av = fl 2 q av where fl 2 is the conformal factor and is the metric of 
asymptotically flat stationary black hole [15l [Itl [17] . A natural generalisation of the surface gravity kdh can be 
introduced by the following relation at the conformal Killing horizon [15l Il6j , 



(3.20) 
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where £ is a conformal Killing vector satisfying C^g^ v — In SI 2 ) g^ v , and the conformal Killing horizon is the 
hypersurface where £ becomes null. Under the conditions Q — > 1 (or constant) and £^£„ — > — 1 at null infinity, Sultana 
and Dyer conjectured that the temperature of conformal stationary black holes Tsd is given by [IB] 

Tsd := (kdh - C t In f! 2 ) . (3.21) 
Tsd is constant everywhere on the conformal Killing horizon, while kjjh is not. On the other hand, Jacobson and 



Kang independently introduced a generalised surface gravity kjk as [17 



V M (C£ a ) = -2 K j K ^. (3.22) 

kjk is invariant under the conformal transformation, while kdh is not. Then under the conditions Q — > 1 and 
£ M £„ — > — 1 at null infinity, they conjectured that the temperature of conformal stationary black holes Tjk is given by 

Tjk := ^ ■ (3.23) 

In fact, it can be shown that the relation Kjk = kdh — C^lnfl 2 holds [17| . and hence Tsd = Ttk =: Tjksd- Thus, 
although Sultana-Dyer [l6| and Jacobson-Kang [T3| considered independently the surface gravity and temperature of 
conformal stationary black holes, they reached the same conjecture. 

For the Sultana-Dyer black hole, the conjectured temperature Tjksd has already been calculated [T3] ■ The conformal 
Killing vector is £ = d n and the conformal Killing horizon coincides with the event horizon r = 2 M. The norm of £ is 
£ M £„, — > —a 2 7^ — 1 at null infinity, not satisfying the unit norm condition for Tjksd- Sultana and Dyer still assumed 
in [12| the temperature Tjksd should be assigned to the Sultana-Dyer black hole. Then, substituting £ into eq. (|3.20[) 
or (f3~2"2"]) , they obtained k jk = k dh -Ajr\ = I /{AM) and 

T JK sd = ^, (3.24) 

which is equal to the Hawking temperature of the Schwarzschild black hole of mass M, but not to our effective 
temperature T c s in eq. (|3. 19|) . However, the temperature of black hole should be given based on the spectrum and/or 
the power of the Hawking radiation. Here, we propose that the physically reasonable temperature of black holes which 
are conformal stationary and asymptotically dynamical is not Tjksd but 

T eS := . (3.25) 

So the effective temperature depends on space and time through the conformal factor. This might be understood in 
an analogy with Tolman's law for thermal equilibrium in the presence of a gravitational field [181 ]. 



IV. SUMMARY AND DISCUSSIONS 



We have calculated the intrinsic power of the Hawking radiation from cosmological black holes for two cases, no 
mass accretion and significant mass accretion. 

For no mass accretion case, we have considered the Einstein- Straus black hole. Our result Tes(4D) in ec l- (|2.49p 
indicates Tes(4D) < Ph(4D)j i-e., the black hole evaporation is suppressed by the cosmological expansion. The ratio 
-Pes(4D)/-Ph(4D) (< 1) is given in terms of e 1 / 3 where e is the ratio in size of the black hole to the cosmological horizon. 
The first correction term is 0(e 1 / 3 ) and therefore currently as small as lO _5 (M/lO 6 M0) 1 / 3 (i/14 Gyr) -1 / 3 , but could 
be significant for the formation epoch of primordial black holes. The evaporation time is essentially the same as that 
of the Schwarzschild black hole as long as its mass is greater than the Planck mass. Furthermore, by comparing the 
functional form of Tes(4D) with that of Ph(4D) in eq. (|2.47p . we can see that the Einstein-Straus black hole has not 
settled down to thermal equilibrium in a finite cosmological time. 

For the significant mass accretion case, we have considered the Sultana-Dyer black hole. This has very different 
properties. Our result Psd in eq. (|3.18[) indicates that the Sultana-Dyer black hole does not evaporate away. Further- 
more the Sultana-Dyer black hole can be regarded as an object in quasi-equilibrium, since the effective temperature 
T ff can be assigned as eq. (I3.19|) . The intrinsic power Psd of the Hawking radiation is consistent with the Stefan- 
Boltzmann law for a black body with temperature T e g. We propose a new definition (|3.25|) for the temperature for 
general conformal stationary black holes. 
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Finally we try to interpret Pes(4D) for the Einstein-Straus black hole in an analogy with quantum radiation of a 
slowly moving mirror in a flat spacetime. For simplicity, we consider a moving mirror x — x(t) in the two dimensional 
Minkowski spacetime with the Cartesian coordinates (t,x). Then it is well known that the moving mirror emits 
quantum radiation of a massless scalar field <f>. When an observer at rest is in the region x > x(t), the observed power 
P m i T of quantum radiation from the mirror is given by [H, [lj| 



p ._ (rn{x)\ 1 Vl-V 2 tfamir . . 

mlr -~ [ W )mil - 12tt (1 - vf dt let ' 1 ' 

where (T^') m j r is a tetrad component, t rct is the retarded time when the observed particle of <j> was emitted from the 
mirror, v — dx(t) / dt\ trot :— x(t rct ) is positive when the mirror is approaching towards the rest observer, and a m i r is 
the proper acceleration of the mirror. In the Minkowski spacetime, the observed power P m i r is equal to the intrinsic 
power of quantum radiation from the mirror. Here we consider the case that the mirror moves slowly, i.e., \v\ <C 1. 
Then the power P m i r is given as 



12tt 



(l + 2v) + 0(v 4 ). (4.2) 



Note that the kinematic effect comes in the radiation power in the form of (1 + 2v) in the lowest order. On the other 
hand for the Einstein-Straus black hole, the relative "velocity" ves of the junction surface E to the black hole at the 
retarded time may be written as 

ves = —rs a ret , (4.3) 

where a ro t := a(t re t)) arid the minus sign means the increase of the relative distance. Equation (|2.9b[) with k = 
gives d 3 et = (2MH Ie t)/r^ — e/rf,, i.e., ves — — e 1 ^ 3 - Therefore, the kinematic correction factor in eq. (|4.2j) coincides 
with the cosmological correction factor in eq. (I2.49|) up to this order. Hence, we can interpret the correction factor 
(1 — 2 6 1 / 3 ) in Pes(4D) is some kinematic effect from the cosmological expansion in an analogy with radiation from a 
moving mirror. 
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APPENDIX: VACUUM EXPECTATION VALUE OF QUANTUM STRESS-ENERGY TENSOR 

1. Two dimensional case 



It has already been recognised for a few decades that many different methods of renormalisation give equivalent 
results (see for example, chapters 6 and 7 in Q). We consider a minimally coupled massless scalar field 0, whose 
stress-energy tensor is given by 

T^u = (j>,fj,<t>,u - -^g^u4>,a4> ,a ■ (A.l) 

The background spacetime is described in double null coordinates (w, v) as 

ds 2 = -D(u,v)dudv . (A.2) 

The field (f> satisfies the Klein-Gordon equation = 0. When a coordinate system (not necessarily null) is specified 
to describe the differential operator □, we can find a complete orthonormal set {/^} for arbitrary solutions of □<£ = 0, 
where u> denotes the frequency of the mode function. The positive frequency mode is the mode function f u which is 
constructed to satisfy the conditions, lu > , (f u , f u i) = S(uj — ui') , (f w , /*,) = and (/*, /*,) = —S(u> — a/), where 
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(/, g) is the inner-product defined from the Noether charge of time translation of <p and /* is complex conjugate to f u , 
called the negative frequency mode. In two dimensional spacetimes, the positive frequency modes can be decomposed 
with respect to the direction of propagation. In the double null coordinates, they are f u {u) = exp(— i ou u) / \/ 4irui 
-iujv)/y/imj. 



and fu(v) = cxp(— i lu v ) / \/ Attlu . Then 
positive and negative frequency modes as 



In the double null coordinates, they are f u {u) = exp( 
the quantum operator cj> is expanded by the complete orthonormal set of the 



du [ a u f u (u) + a t /* („) + f u (v) + fet /•(„)" 



(A.3) 



The canonical quantisation presumes the simultaneous commutation relation between <fi and its conjugate momentum, 
so that {flu} and {b^} are harmonic operators satisfying the commutation relations; [a u , a^,] = 6(lu—lu') and b^,] — 
5(w — lu') and all others vanish. They define the Fock space of quantum states and give particle interpretation. The 
vacuum state |vac) is defined as a quantum state satisfying a w |vac) = 6 w |vac) = for all lu. 

If we choose different coordinates (u,v), a natural orthonormal set of mode functions is {fa}, where fuj(u) = 
exp ( — i lu u)/y/4iru> and f u (v) = exp(— i lu v)/\/4mjj. Then the expansion (| A.3|) gives different harmonic operators 
{a w } and {b u }. These harmonic operators define another vacuum state |vac) |vac) ) if there arises the mixing of 
positive and negative frequency modes (f u , f*,j ^ between the two coordinate systems. Thus, even if a quantum 
state is initially set to be a vacuum state, this does not remain vacuum but corresponds to an excited state associated 
with the coordinate system natural to an observer at the final time if the mixing of positive and negative modes arises. 
This will be interpreted as quantum particle creation in curved spacetimes. 

The quantum expectation value of the stress-energy tensor (vac|7)„,|vac) is calculated by substituting the quantum 
operator (|A.3|) (after replacing a u and b w with and b u ) into the stress-energy tensor (|A.1|) . However, (vac|T M ^|vac) 
diverges even for flat background cases. Therefore, we need to renormalise the stress-energy tensor. We do not get 
into the details of the regularisation method but only quote the result [20j , 



K 

48^ 



9\iv 



(A.4) 



where (Tp P ) is the renormalised expectation value of (vaS|r^p|vac), 1Z is the Ricci scalar of the background spacetime, 
and 0p P is a symmetric tensor whose components in the coordinate system (u, v) on which the vacuum |vac) is defined 
is given by 



1 
1 

= Ovu = 



3 

2\D 



D 



D 

. VV 

D 



(A.5a) 

(A.5b) 
(A.5c) 



where D(u,v) = —2guv- The renormalised expectation value (T^ v ) of (vaclT^lvac) in the other coordinates (u,v) is 
calculated from the above components through the usual coordinate transformation for tensor components, 



3x^ dx v 



(A.6) 



2. Four dimensional case 



The renormalised vacuum expectation value of stress-energy tensor (T M „) of some matter field in four dimensions 
may also be calculated along with the canonical quantisation formalism as shown for the two dimensional case in 
previous section. However the path integral quantisation formalism is more convenient to summarise (T^) on a four 
dimensional conformal spacetime. The effective action W of a quantum matter field cf> on a spacetime of metric g^ v 
gives the vacuum expectation value of quantum stress-energy tensor. W can be evaluated by the path integral method 
and the vacuum state |vac) is specified by the Green function of (f> used in evaluating the path integral. However 
the precise path integral form of W is not important here. W is decomposed into two parts as W = W rcn + Wdiv, 
where W Tcn is the renormalised part and Wdiv is the divergent part. The functional differentiation of W Ien gives the 
renormalised vacuum expectation value (T^), 



(A.7) 
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We consider the case that the metric g^ v is conformal to the other one as 

g^ = n 2 g^ u , (A.8) 

and the matter field is a conformally coupled massless scalar field satisfying (□ — lZ/6)<fi = 0. On the other hand, 
we get by definition of functional differentiation, 

Wrcn - W Tcn = J £pS Sg^ d*x , (A.9) 

where W ren is the renormalised effective action obtained from W Ten with replacing g^ v by g^- Then considering the 
functional differentiation only by the conformal transformation, 5g^ v = —2g f * L/ O -1 SCI, the effective action is expressed 
as 

W Icn = W Icn - / g a0 {T aP ) —^j^x . (A. 10) 

Substituting this into eq. (IA.7j) . we get 

(T Al ,) = -^(f^)--L-^ ; [g"P{T a p) S -£-V=9d*x, (A.ll) 



W ^' V^9Sg^ J y x ap/ Q 
where 6^ = g^ a g av , g^ag aa — and the general relation, 

Sg au bg av 

are used to get the first term of the right-hand side of eq. (|A.11[) . The trace g a ^{T a p) is usually called the conformal 
anomaly or the trace anomaly, and it is well known that the divergent part Wdw gives the conformal anomaly as (see 
§6.3 in [6j for example) 

g (T a0 )- — -^-. (A.13) 

Hence substituting this expression of the conformal anomaly into eq. (|A.11[) and using eq. (|A.12[) and eq. (|A.9[) with 
replacing W xen by Wdw, we obtain 

The divergent part Wdiv can be evaluated from the Green function of the matter field <j>. We do not follow the 
details of the calculation of Wdiv> but quote only the result for (Tu V ) for the conformally coupled massless scalar field 
<f> on the spacetime of metric (|A.8|) (see [2l[ or §6.2 and §6.3 in for detail), 



(V) - ^2 ( T ^> 2g80 ^ 2 ( 6 Y^j + 2g80 ^ 2 (qX^ Y»u ) , (A.15) 



where 



X„ v := 2 V l S7 v n~2g^un+^K 2 g^-2TZR^, (A.16a) 

Ypv := -R£R av + \llR llv + ^R a pR a Pg^-^n 2 g liV , (A. 16b) 

where i?^„ and 1Z are the Ricci tensor and scalar with respect to g^ respectively, and and are defined 
similarly with respect to the metric g^. Equation (|A.15j) is the generalisation of eq. (6.141) in Q to the general 
conformal spacetimes of metric (|A.8[) . 
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